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Abstract 

The Hodrick-Prescott filter represents one of the most popular method for trend- 
cycle extraction in macroeconomic time series. In this paper we provide a multivariate 
generalization of the Hodrick-Prescott filter, based on the seemingly unrelated time 
series approach. We first derive closed-form expressions linking the signal-noise matrix 
ratio to the parameters of the VARMA representation of the model. We then show that 
the parameters can be estimated using a recently introduced method, called “Moment 
Estimation Through Aggregation (META)”. This method replaces the traditional 
multivariate likelihood estimation with a procedure that requires estimating univari¬ 
ate processes only. This makes the estimation simpler, faster and better-behaved 
numerically. We prove that our estimation method is consistent and asymptotically 
normal distributed for the proposed framework. Finally, we present an empirical ap¬ 
plication focusing on the industrial production of several European countries. 


1 Introduction 


The extraction of trend and cycle components from economic time series represents an 
important tool for economic analysis. Several univariate methods have been discussed 
by the literature (see for example Beveridge & Nelson (1981), Watson(1986), Harvey & 
Jaeger (1993), Canova(1998), Baxter & King(1999), Harvey & Trimbur(2003)). One of the 
most popular method, widely employed in macroeconomics, is the “smooth-trend model”, 
generally known as Hodrick-Prescott(1997) hlter. As a matter of fact, this method was 
suggested a long time before by Leser(1961) for trend extraction (see the discussion in 
Mills(2009)). This approach extracts a stochastic trend which moves smoothly over time 
and is uncorrelated with the random irregular term representing the cyclical component. 
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The ratio between the variances of the two noises (i.e. the signal-to-noise ratio) is the key 
scalar that determines the “smoothness” of the extracted trend. For example, Hodrick &: 
Prescott(1997) suggest specific values for time series observed at different frequencies. 

Despite the recognized importance of this method in empirical analysis, we still know 
little about the multivariate case. This represents a relevant framework since it allows 
extracting multiple trends that might share similar dynamic behaviors (such as the case of 
common trends as in Stock & Watson (1988)). A notable exception is Kozicki(1999) who 
discusses the multivariate extension of the Hodrick-Prescott filter. However, the author 
assumes the same single common trend for the whole system of equations. This assumption 
is rather restrictive and is relaxed here. 

We provide analytical results for the Hodrick-Prescott filter in the multivariate case. 
More specifically, we derive closed-form results linking the signal-noise matrix to the pa¬ 
rameters of the VARMA model representing the stationary representation of the structural 
model. Similar expressions for some trend-cycle models can be found in Morley et al.(2003) 
and Oh et al. (2008); however, none of the mentioned papers deal with the multivariate 
case. Establishing explicit relations between the two forms is relevant since one can de¬ 
rive/estimate the structural parameters from the reduced form ones (and viceversa). 

Relying on these relations, we build a fast and simple estimation method for the 
covariance matrices and extract from it a change of variable matrix that decouples the 
model into d uncorrelated ones with the same form. These models can then be estimated, 
each with its own optimal signal-noise ratio. Our method generalizes the so-called META 
approach of Poloni & Sbrana(2014), which was initially developed for the multivariate 
exponential weighted moving average model. We prove in the appendix that the resulting 
estimator is consistent and asymptotically normally distributed. 
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Finaly, we show an example for the practical use of our results in extracting the trends 
from the industrial production series of some European countries. 

2 Theoretical results 

The unobserved components representation of the smooth-trend model, also known as 
Hodrick-Prescott filter, was firstly used by Akaike(1980). Here we consider the multivariate 
state-space representation (for the univariate case see Harvery & Trimbur(2008)) 

yt = + et, 

Ht+i = l^t + /3t, (1) 

Pt+i = Pt+it, 

where the vectors yt, yt, it are of dimension d. In addition t = 1,2, ■ ■ ■ ,N represents 
the number of observations. Contrary to Harvey & Trimbur(2008), here we relax the 
assumption of normality of the noises. We assume that the noises et and it are i.i.d. with 
zero mean and covariances 


e* 


1- 

o 

w" 

A 


0 


where and 0 are d x d matrices. In the scalar case the ratio between the variances 

of the two noises (i.e. the signal-noise ratio) is the key scalar that determines the 

“smoothness” of the extracted trend ft. For example, the smaller the ratio, the smoother 
the extracted trend (see Kaiser & Maravall (2001)). In the multivariate case the ratio is 
a non-symmetric matrix (in general) whose eigenvalues play the role of the scalars that 
determine the smoothness of the trends of the transformed series (this will be clarified 
below). 

The nonstationary state-space system is called a structural process and the co- 
variances in m are called structural parameters. Its stationary representation zt = 
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yt — 2yt-i + yt -2 is a (second order) integrated vector moving average of order two (See 
Harvey(1989) and Maravall & Del-Rio(2007) for the univariate case). Using Wold repre¬ 
sentation theorem, we can reparametrize it as 

Zt = = (/ — 2L -|- L?‘)yt = Ct -2 {I ~ -|- L‘^)et = 


( 3 ) 


= {I + QiL 02L^)r/t, with E{r]t'n't) = 

where L is the backshift operator and I is the d x d identity matrix. This form is known 
as reduced form, with parameters 0i, 02 and U. 

The autocovariances of zt can be expressed as functions of both structural and reduced 
form parameters as follows; here and in the following, we use the symbol M' to denote 
the (conjugate) transpose of a matrix M. 


To = E {ztz't) = 6S, + Eg = U + 0iU0; + 02U0'2 (4a) 

Ti = E {ztz't_^) = -4E, = 0iU + 02^0;, (4b) 

r 2 = ^ {ztzi_^) = E, = 02U, (4c) 

Tj = E {ztzt_j) = 0, for |jj > 3. (4d) 

Note that the autocovariance matrices are all symmetric i.e. r _2 = T^ = r 2 and 


r_i = = Ti; this is a characteristic feature of this model that we shall exploit in 

our computations. 

Computing 0i, 02 and Q from the covariances T^ in (j4|) (or, equivalently, from E^ and 
E^) requires solving a system of nonlinear equations. We wish to show how the solution 
can be determined explicitly in closed form. 

One can gather together the autocovariances to form the autocovariance generating 
function (ACGF) (see Harvey(1989)), which is a rational function in the formal variable 
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L 

r(L) := + TiL + Tq + 

= - 4S,L + 6S, + S. - 4S,L-^ + + E.. 

The relation among the autocovariances and the MA(2) parameters can be embedded 
in the following factorization 

r(L) = (02^2 + 0iL + i)n{i + e[L-^ + q' 2 L-^). (5) 

If the matrix polynomial 02-^^ + 0iL + I has no roots inside the unit circle, the MA 
process is called invertible. If the ACGF is such that r(z) is invertible for each z on the 
unit circle {z G C: \z\ = 1}, such factorization (called eanonical faetorization) exists and 
is unique; for a formal proof of this statement in the matrix case, see Gohberg, Lancaster 
and Rodman(1982), Theorem 4.1 (and the following remark). Therefore, if we determine a 
canonical factorization ([5]) of the AGGF, then its coefficients 0i, 02 and Q, must coincide 
with the MA(2) parameters. 

0ur strategy is constructing explicitly such a factorization of the AGGF. We start from 
the scalar case, then move on to the multivariate one. 

2.1 The univariate case 

Here we provide the algebraic linkage between the scalar signal-noise ratio and the moving 
average parameters. This linkage is well-known and is discussed for example in McEl- 
roy(2008), for the HP-filter, as well as in Sbrana(2011) for the generic local linear trend. 
The result in Proposition [1] is instrumental for the multivariate case. 

Proposition 1. Let zt follow an univariate version of the model 0, and suppose > 0, 
Eg > 0. Let 6 = E^EF^. Then, zt follows an invertible MA(2) zt = {1 + OiL + 62 L?‘)ut 
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process with coefficients 


2 4 + 6/1 4 + V-26 + 2 V (^2 + 16^ 


and oj = Var ut = @2 

Proof. In the scalar case, the autocovariance generating function takes the form 


7(L) = 


(L - 1)^ 


^2 

The complex zeros of 7 (T) can be determined through the following process 
(s-1)^ 


(s-1)^ 


= --5, 

= ±iVd, 


0 = s"^ - (2 ± i\/dj s + 1 
2 + it ± V4it — 


s = 


, where t = ±\/5, 


( 6 ) 

(7) 

( 8 ) 


where the last step is simply the formula for the solution of a quadratic equation. The 
possible choices of the ± signs in ([ 8 ]) give the four solutions. Each of the two quadratic 
equations in ([T]) has two solutions with product 1 and sum 2 iLi\/S, by the roots-coefficients 
relations. Since their product is 1, one of them lies inside the unit circle and one lies outside 
(they cannot have both modulus 1 , otherwise their sum would have modulus at most 2 , 
which is in contradiction with |2 ± i^/ 6 \ > 2). Moreover, since Q is a real equation, we 
know that the solutions come in conjugate pairs. Putting all together, we have proved 
that the solutions returned by the formula ([5]) can be written as (s, s, 1 /s, 1 /s) for some 
complex number s with |s| < 1 . 

To respect invertibility, the polynomial 1 + 61 L + 02 must have as its roots the two 
roots with modulus larger than 1, i.e., 1/s and 1/s. Hence, 

1 + 0iL + 027 ^^ = (1 - (s + s)L + ssL^). 
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To avoid troubles with the signs, we derive an equation for 6i directly. The univariate 
equivalent of ([5]) is 

7(T) = {02L^ + 6iL + l)w(l + 9\L ^ + 92 L ^). (9) 

Equating coefficients in Q, we get the following system of equations 

CJ, = 72 = 920J, 

-4cj, = 7 i = 6 'i( 6'2 + 1 ) 0 ;, (10) 

6cj£ + = 7o = ^(1 + + ^ 2 ) 

(since we are in the scalar case, for the sake of clarity we replaced each uppercase letter 
with the corresponding lowercase one). We can easily derive 

= uj = 92^a^, ( 11 ) 

and use these relations to eliminate variables and get a single equation in 9i. After some 
computations, we obtain 

9\ + S9l + (24 + d)9l + (32 + 46)9i + 16 = 0, 

whose four solutions are given by choosing ± signs in 

-2 ± -26 ±2^6^ + 166. 

Only one of these solutions corresponds to s + s (the other ones being ^ + + ^ + 

which would give rise alternative non-invertible MA(2) representations). It is easy to tell 
which one is correct: it should be real, and this implies that we have to take the + sign 
on the right; and since |s + s| < 2 we need the plus sign on the left as well. 

Finally, one can use (fTT]l to get back 9i and w. □ 
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2.2 The multivariate case 


As claimed in the introduction, it is possible to make a linear change of variable that 
transforms ([T]) into d separate uncorrelated processes. We show its form explicitly in the 
next result. 


Proposition 2. Let zt follow the model and suppose that and are positive 
definite. Let = M'M be a Cholesky decomposition, and = QAQ', with 

QQ' = I and A = diag((ii, (^2, ..., <5^); he an eigendecomposition, and let P = M'Q. 

Let moreover 


1 I / -aj. 4 - 

Ofc = + —\ —26k + + 16(5fc, fik = -=- ; 

2V ^ 4 + afc / r—- 


4 + —26 k + + 16(5fc 

Then, the unique invertible VMA(2) representation of zt is given by 


01 = Pdiag(ai,a 2 ,...,ad)P \ 02 = Pdiag(/3i,/32,...,\ (12) 

n = Varut = Pdiag(/3]"\/32'\ ... = 0^^S,. (13) 


Remark 1. When Eg is positive semi definite with some roots equal to zero then the model 
contains common trends. This is the case when cointegration arises. In this case some of 
the 6k are equal to zero and therefore the corresponding ak and fik equal to -2 and 1 
respectively. 

Notice that EgS^i = M'(M')"^SgM-i(M')"^ = M'QAQ'{M')-^ = hence 

P is an eigenvector basis of SgE^^. 0ne may wish to choose P as an arbitrary eigenvector 
basis of SgS^^ instead of the more complicated definition in the theorem. In that case, 
then (I12p and Ll = still hold, while the other equality in (jl3D may fail. 










Proof. The autocovariance matrices of zt in (l3|) are 


r2 = E {ztz'i-_2) = Se, 

Ti = £' {ztz'i._i) = -4S,, 

Tq = E [ztz[) = 

In addition, r _2 = T^ = r 2 and r_i = = Ti, and all the other autocovariances Tj with 

|i| > 2 are zero. 

Since only the two symmetric matrices appear in these expressions, one can hnd 

a change of variables that decouples the components of zt- Namely, we set zt := P~^zt, in 
this way, it is easy to verify that P~^E^{P')~^ = I and P~^E^{P')~^ = A, and thus 

t2 = E {ztz[_2) = I, 
ri = E {zt~z't-i) = -4/, 
ro = E {ztz[) = 6/ + A, 

and again f _2 = T '2 = P 2 and f_i = = fi, while ffc = 0 for |A;| > 2. 

Hence each of the components of zt follows a scalar MA(2) uncorrelated from the other 
components, with = 1 and = 5^, where k is the index of the component. Using 
Proposition [H one can thus derive 

{zt)k = (1 + akL + (5kL^){ut)k, Var(ht)fc = 

Finally, we undo the change of variables used to dehne zt, to obtain that zt follows the 
VMA(2) process (fT^ - (fT^ . where ut := Put- □ 

Therefore it is possible to reconstruct the parameters of the reduced form of the process, 
in closed form, by knowing only its autocovariances. In empirical analysis, one might 
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be tempted to use the sample covariance estimates such as fg = N~^ and 

f 2 = N~^ ^tz't _2 but this is generally not recommended due to lack of accuracy. 

In the next section, we provide a method that allows estimating more accurately Eg = 
r 2 and Sj = To — 6 r 2 . 

3 Moment estimation through aggregation (META) 

The closed-form results obtained in the previous section are relevant in empirical analysis 
if we have accurate estimates of the autocovariances of ([3]). Indeed, these allow one to 
reconstruct the signal-noise matrix ratio and therefore extract the multiple trends from 
system m- In this section, we describe an estimator of the autocovariances that general¬ 
izes the algorithm of Poloni &: Sbrana(2014) to this model. The main cost of the estimation 
procedure consists in estimating several univariate constrained MA(2) models of the form 
(ttH); this makes it quick and practical even in cases of large dimension, without conver¬ 
gence issues. We then combine these estimates to obtain the desired autocovariances T^. 
Beside these practical advantages, the reader should be aware that the META approach 
does not guarantee to yield positive definite and This is especially true in small 
samples and also when we approach the cointegration case (that is when the one or more 
roots of are closed to zero). This issue and a simple proposal to fix it are discussed 
below. 

The estimation procedure can be described as follows. Let 

W := {cj: 1 < i < d} L) {ei + ej: I < i < j < d}, 
where denotes the Tth column of the identity matrix 

1. For each of the vectors u) G W, construct the aggregated scalar process 
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:= w'zt, and estimate it using a maximum likelihood estimator. 

2. Using the estimated parameters of the process compute autocovariances 
{not sample autocovariances!) and construct the (scalar) ACGF 7^ ^\L) of x^^\ 


3. 


4. 


As the ACGF F(L) of yt is symmetric, it is possible to reconstruct it in closed form 
by knowing = w'T{L)w for each w E W. Indeed, by Poloni & Sbrana(2014), 

Lemma 5, 




(14) 


Use the method in Section 2 to compute the signal-noise matrix ratio as well as the 
parameters of the VMA representation of zt- 


As can be seen easily from the formula 7 ^'^^ (L) = w'T{L)w, the aggregated scalar processes 
are MA(2) with = — 47 ^*^^ As in (lll|) . we can obtain from this relation 


q(w) _ "1 


4 + 0) 


(w) ■ 


(15) 


We estimate them using a constrained maximum-likelihood estimator enforcing ()15p . The 


procedure produces 0^ , 02 = —and . We can then use the formulas (cfr. (flUIl l 

72 — ^2 7 


7 , = 05 “') ( 0 ("') + l)cD(-) = -472, (16) 

7 o = J;(-)(1 + ( 0(-))2 + ( 0 (-)) 2 ). 

Notice that the estimated autocovariances F^. of the multivariate process will satisfy au¬ 
tomatically the constraint Fi = Fl^ and Fi = —4F2, by the linearity of (1141) . If one 
were to use a direct multivariate maximum-likelihood estimator, these conditions would 
be harder to represent explicitly as a restriction on the parameters 0i and © 2 , and hence 
troublesome to estimate. 
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Finally, it can be shown that the META estimator is consistent and asymptotically 


normal distributed. Indeed, the following theorem holds. 


Theorem 1. Consider the META estimator described in Section\^for the reduced pa¬ 
rameters 01 and 02 of the invertible, stationary and ergodic process © with i.i.d. noises 
with variances ©• The estimator is consistent. If, in addition, the fourth moments of 


the noise vector 


are finite, then the estimator is asymptotically normal. 


A proof can be found in the Appendix. It is relevant to remark the advantages of using 
the META estimator. Indeed the standard maximum likelihood estimation of the system 
is not trivial especially for medium-high dimensional systems. Similar conclusions 
hold for the likelihood estimation of the reduced form ([21). By making use of a univariate 
estimation approach, the META provides a simple estimation alternative, based on the 
likelihood principle (providing accurate estimates), much faster than the full multivariate 
likelihood approach. The procedure can be implemented by standard packages since it 
does not adopt any sophisticated maximization algorithm. 

As already mentioned above, the major issue with this numerical procedure is that 
there is no guarantee that the estimated process yields positive definite values of the 
structural parameters = f 2 and = Tq — 6 f 2 . This is especially true for Eg, since it 
is derived as the difference of two estimated matrices and hence it might suffer from error 
accumulation. 

A simple fix for this issue is enforcing positivity by adding a suitable multiple of the 
identity to E^ and Eg when necessary. Ultimately, this can be justified by the assumption 
that our data come from a model of the form ®, thus their positivity is a modelling 
requirement. We describe this regularization procedure in more detail for our practical 
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example. 


3.1 Trend extraction 

The key factor for trend extraction using the Hodrick-Prescott filter is the signal-noise 
ratio. In the scalar case, the ratio can either be chosen arbitrarily (as proposed in Hodrick- 
Prescott (1997)), or can be estimated using the result in Proposition [TJ This choice is also 
present in the multivariate case. In the vector case, we have shown how to estimate the 
signal-noise matrix ratio = ToTg ^ — 6/ = PAP~^. Each of the diagonal elements of 

A represents a scalar signal-noise ratio that can be used to hlter the associated component 
of the transformed system yt := P~^yt, whose VARMA representation is diagonal. One 
can then obtain the smooth-trends of the system yt by premultiplying the transformed 
trends yt by the matrix P. 

A possible variant is choosing arbitrarily the values of the signal-noise ratios on the 
uncorrelated processes yt rather than on yt. This produces a trend extraction method 
that, while still choosing SNRs arbitrarily, keeps into account the fact that the covariance 
matrices ([2]) are in general not diagonal and the processes are correlated. 

A direct multivariate approach for trend extraction, that does not diagonalize the 
system, is suggested by McElroy &: Trimbur(2015). These authors provide also results for 
the trend model as in ([T|) when has reduced rank. 

4 Detrending multiple time series: an empirical application 

Here we provide an empirical example dealing with the industrial production of some ma¬ 
jor European economies. The time series data, relative to the total industrial production 
by country, have been downloaded from the online statistical database of the Organisa- 
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tion for Economic Co-operation and Development (2014). We employed monthly season¬ 
ally adjusted observations (expressed as index numbers) for the following eight countries: 
Belgium, France, Germany, Italy, Netherlands, Portugal, Spain, United Kingdom. The 
sample runs from May 1974 until March 2014 for a total of 479 observations. 

As described in Section [3l the bulk of the estimation cost consists in estimating 
= 36 scalar constrained MA(2) models, followed by some quick computations. For 
our experiments we chose to use Wolfram Mathematica v9, since that software contains a 
ML estimator general enough to allow enforcing the constraint (|15l) explicitly. Note that 
the scalar estimator used is unconditional maximum likelihood. 

The estimated covariance matrix Sg had a small negative eigenvalue ~ —1.5 x 10“^. 
Indeed, as we noted above, there is no guarantee that the method produces positive 
definite estimates in empirical examples. To address this issue, we regularized the estimate 
by adding a suitable multiple of the identity to i.e., -|- al^. We chose 

a = 0.0015428533. This number is chosen so that the smallest eigenvalue of equals 

= 0.0000694, which is the standard signal-noise ratio recommended in software 
packages such as EViews for Hodrick-Prescott filtering of monthly data. 

We obtain the following estimates for the covariance matrices of the noises ([2]) 


0.8429 

0.02121 

0.7632 






\ 

0.09744 

0.1419 

0.7497 






0.1871 

0.06468 

0.06268 

1.159 





-0.1544 

0.2812 

0.06318 

0.1073 

2.583 




0.06602 

0.2387 

0.2089 

0.265 

0.1438 

3.048 



0.06847 

0.247 

0.08804 

0.131 

0.1751 

0.3083 

1.686 


0.09582 

0.1087 

0.0776 

0.1215 

0.1535 

0.1928 

0.1191 

0.6269 / 
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0.07748 

0.06312 

0.05871 






\ 

0.08688 

0.08175 

0.1239 






0.0734 

0.06911 

0.09324 

0.1618 





0.04262 

0.03635 

0.05105 

0.05275 

0.03846 




0.02446 

0.01703 

0.02533 

0.03828 

0.01042 

0.02073 



0.04151 

0.04143 

0.05226 

0.06523 

0.02957 

0.01757 

0.04876 


0.03398 

0.03357 

0.05591 

0.03005 

0.02341 

0.005987 

0.01746 

0.04468 / 


The META estimates are very close to those produced by STAMP 8.2 (Koopman, 
Harvey, Doornik &: Shephard, 2007), which employs a maximum likelihood estimator for 
the structural system. Indeed, a comparison between the META and the STAMP output 


shows that the relative errors in the Frobenius norm (root mean squared error of the 


matrix entries) are 


^META 


_ y]STAMP|j 


Y^META _ ySTAMP 


YSTAMP 


^^META _ QSTy4MP|| 
||H^TAMP||^ 


0.073, 

0.177, 


0.074, 


(17) 

(18) 


(19) 


where H-Hp. is the Frobenius norm. 

However, while the estimation in STAMP takes about 40 seconds, the one with META 
takes just 14 seconds. This shows the clear computational advantages of using our simple 
estimator. Using the estimated matrices, the signal-noise matrix ratio is derived using the 
procedure described in Proposition [2j We have 


■ 0.4601 

0.08478 

-0.4455 

0.4898 

-0.063 

-0.2311 

0.109 

0.3412 

0.4225 

0.07566 

-0.1639 

-0.1915 

-0.134 

0.05911 

0.3999 

-0.578 

0.6077 

0.2573 

0.0001279 

-0.3538 

0.2624 

-0.0273 

-0.342 

0.04939 

0.5831 

-0.7558 

0.4043 

0.2669 

0.05411 

-0.07511 

0.06081 

0.03422 

0.2834 

-0.01842 

-0.03882 

0.4361 

-0.543 

0.9852 

-0.7428 

-0.7025 

0.1505 

-0.1733 

-0.09401 

0.2449 

0.11 

-1.192 

-0.6928 

-1.007 

0.3025 

-0.2091 

-0.08031 

-0.4969 

-1.071 

-0.3526 

-0.1477 

0.0638 

0.259 

0.3946 

0.4872 

0.2695 

-0.2398 

-0.1559 

0.0672 

-0.08842 
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and 


/ 0.3127 0. 0. 0. 0. 0. 0. 0. 

0. 0.08487 0. 0. 0. 0. 0. 0. 

0. 0. 0.03726 0. 0. 0. 0. 0. 

. _ 0. 0. 0. 0.01205 0. 0. 0. 0. 

^ “ 0. 0. 0. 0. 0.01085 0. 0. 0. 

0. 0. 0. 0. 0. 0.0054 0. 0. 

0. 0. 0. 0. 0. 0. 0.004519 0. 

V 0 . 0. 0. 0. 0. 0. 0. 0.00006944 

The values of A are then used to extract separately the trend from each series of the 
system yt ■= P~^yt- This univariate procedure can be easily carried out with any statistical 
software. Here we make use of EViews version 8. Once these trends are extracted, the 
final step is pre-multiplying these trends by the matrix P in order to obtain trends for the 
original multivariate system yt- 

Here we compare our empirical results obtained using the META approach with two 
standard univariate approaches. The first approach is the ARIMA-model-based (AMB) 
approach as suggested by Kaiser & Maravall(2005). The AMB approach can be considered 
as the univariate analogous of the META approach since the signal-noise ratio is estimated 
(rather than fixed apriori) by employing an ARIMA(0,2,2) model separately for each series 
of the system. This procedure is implemented in Eviews through the function SEATS 
(Signal Extraction in ARIMA Time Series) by performing an ARIMA-based decomposition 
of an observed time series into unobserved components. The SEATS algorithm in Eviews 
was developed by Victor Gomez and Agustin Maravall. 

The second approach, the most employed in standard practice, consists of fixing the signal- 
noise ratio to a prescribed value (see for example Hodrick & Prescott (1997)). For monthly 
series EViews suggests a signal-noise ratio of = 0.0000694. The main difference 

between these two univariate approaches is that fixing the signal-noise ratio results in 
extracting a much smoother trend compared to estimating the signal-noise ratio using 
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an ARIMA approach. This in evident in our empirical results. Figures 1-8 report the 
industrial production index (in gray) together with the trends extracted using the META 
approach (thick black line), the AMB approach (dotted black line) and the approach that 
fix the signal-noise ratio (tiny black line). As noted above, the estimation sample for the 
three competing approaches is 1974-2014. However, for the sake of clarity, for each country 
we show the results in two separate charts; one referring to the sample 1974-1994 and the 
other one to the sample 1994-2014. This is done in order to better focus and compare the 
different extracted trends. 

First of all, one can observe that fixing the smoothing constant provide much smoother 
trends compared to the AMB approach. Interestingly, the META approach provides trends 
with mixed level of smoothness. On the one hand, results relative to Spain and UK show 
that the META provides trends that are very close to the AMB approach. This is not 
surprising since both filtering procedure estimate the signal-noise ratio. On the other 
hand, looking at the results for France and Germany, the META extracts trends that are 
very close to the standard practice of fixing the signal-noise ratio (this is especially evident 
in the sample 1994-2014). Finally, the results relative to Belgium, Italy, Netherlands and 
Portugal show that the META approach provides trends that, in terms of smoothness, are 
somehow in between the two univariate approach. Therefore, our empirical results show 
that estimating the signal-noise ratios using the META approach does not necessarily 
deliver different outputs compared with the standard practice of fixing these parameters 
apriori. We believe nevertheless that our procedure is more rigorous since it is based on 
a robust estimation method. On the other hand, imposing a predetermined single signal- 
noise ratio for all series represents a rather simplistic assumption. Indeed, the estimated 
values for the signal-noise ratios are quite far from being constant, in this example. On 
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top of that, there is no general consensus on how to choose the signal-noise ratio (see for 
example the discussion in Ravn & Uhlig (2002)). 

5 Conclusions 

This paper provides closed-form results for the Hodrick-Prescott filter in the multivari¬ 
ate case. In addition, a simple and fast method is suggested to estimate the VARMA 
parameters of the implied structural process. As a consequence, the signal-noise matrix 
ratio can be quickly be estimated and used for filtering the underlying system of equa¬ 
tions. Contrary to the standard maximum likelihood estimation, the main advantage 
of our method is that it is exempt from the numerical and convergence issues of high¬ 
dimensional minimization procedures. Indeed, the META estimation procedure uses only 
univariate model estimations as its computational core. Another significant advantage is 
that these scalar estimations are computationally independent and hence very suitable for 
parallel computation. 

In general, we remark that our results, as well as the estimation procedure, are valid 
when it is assumed that the data generation process ([T|) holds for the entire system of 
equations. Clearly this hypothesis might represent an issue in empirical analysis when 
departing from the assumption of a single dynamics for all the equations. Nevertheless, 
if the aim is simply extracting smooth-trends from multivariate time series, these results 
might simplify considerably the calculations. 
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6 Proof of Theorem 1 


We first recall some definitions. For a sequence A: G Z, if there is a > 0 such that 
||^(*^)|| = 0{p{k)'i(j^) for some polynomial p{k) and for k —)■ oo, we say that g decays with 
rate V’- A causal linear process of an uncorrelated time series {Yt},Yt G is a time 
series of the form 

OO 

gt:=Y.g^’^^Yt_k, 

k=0 

where g^^"^ G Using the backshift operator, one can write alternatively gt = 

and define the so-called transfer function g{L) := If 

decays with rate V’j we say for shortness that the process gt decays with rate f). Often one 

can rewrite a transfer function g{L) as a matrix fraction g{L) = q{L)~^p{L); in this case, 

the process decays with rate max{| 2 ;| : z G C, q{z) = 0}. 

We start by proving a slight variation of Poloni Sz Sbrana(2014), Lemma 10. Here the 

1 / 

notation ||A ||2 denotes the norm of random variables ||A ||2 = IE [A^] 

Lemma 1. Let (li)igz be a sequence of i.i.d. vector-valued random variables in with 
mean 0 and finite fourth moments, and let 

OO 

gf.= J29^''^Yt-k, 

(20) 

ht 

k=0 

be two causal linear functions ofY, with g^^\h^^'i G Let moreover be the a-field 

generated by Yt for a < t < b. Suppose that g^^^ and h^^'i decay with rate ip < 1. Then, 
= gtht — E [gtht] is such that ||.^o ~ E [.^o I II 2 decays (as a function of t) with rate 
Ip, too. 
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Proof. We follow the proof in Poloni &: Sbrana(2014), generalizing it to this case. For a 
fixed t > 0, we may write 


k=0 k>t 


( 21 ) 


■■=Pt 
rO 


:='?t 


where pt is a function in the cr-field and qt is independent from it, and similarly 




( 22 ) 


k=0 


k>t 


\=rt 


:=st 


One has 


E [goho I = E [{pt + qt){rt + St) \ J^_t] 


-- ptvt + E [gt I j] rt+pt'E [st I +E [qtst 

V ^ V ^ 

^ V ^ 


.1 = nxTx -I- 


thus 


||?o - IE [Co I II2 = \\90ho - E [goho \ J^t] II2 = htn + PtSt + qtSt - E [qtSt \ II2 

< Ikt|l2lktll2 + l|Ptll2ll'Sill2 + 2|kt|l2ll'St|l2- 

Since the decompositions ([2T]) . ([22]) are into independent (orthogonal) components, one 
can estimate 


\\Pt\\2 < Il50|l2> 

Iktll2 = 0 1 

Eii^iO 

\k>t / 
/ \ 

= o{if^ + +... 

II 

0 

k’i 2 — ^0 2> 

I|st|l2 = 0 1 

\k>t ) 

= 0{^f^) 



(since Efe>oV’‘^^ = 1 ^ = 

Hence ||Co - E [Co | II 2 = 0(V’*)- □ 
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We are now ready to consider the asymptotic properties of the maximum likelihood 
estimator for the scalar processes We recall that they follow a a scalar MA(2) of the 
form 




(w) I aM (w) 


Jw) 


4 + 0} 




(+v 


= 


(23) 


where is an uncorrelated white noise sequence (but, for wi ^ W 2 , and will 
in general be correlated). 

First of all, we prove that these MA processes admit an invertible representation 
(cfr. Poloni and Sbrana(2014), Lemma 6). 


Lemma 2. Let w G be given, with w ^ 0, and let = w'yt be the aggregate 

proeess of a process yt G with ACGF r(L) such that detr( 2 ;) / 0 for each z on the 
unit circle. Then, the ACGF 7 (L) = w'T{L)w has also no zeros on the unit circle and 
hence x^^'^ has an invertible representation. 


Proof. For z on the unit circle, F( 2 ;) is a Hermitian positive semidefinite matrix (thanks 
to ([5])). So w'T{z)w = 0 can hold if and only if F( 2 ;)t(; = 0, i.e., if detF( 2 ;) = 0, which is 
against our assumption. □ 

Therefore we may safely assume that the polynomials 9{L) = 1 + 6^ L - 

have all their roots outside the unit circle. 

The corresponding negative log-likelihood function is 

N 2 

L{9i,uj) = = ^logui + ^ 


N 


t=l 


where the sequence vt = vt{x^'^'l, 0i) is generated by vt = — OiVt-i + ■^^vt- 2 - As in 

Poloni & Sbrana(2014), we may ignore the issue of the initial data and set vq = u_i = 0 
(quasi-likelihood). Indeed, as in many similar models, the influence of the choice of vq 
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and V-i on the estimator is negligible asymptotically (cfr. Box &: Jenkins(1976), Section 
7.1.3). 

First of all, we check that the exact values of the parameters correspond to 


an isolated maximum of IE 


£['^\9i,(jj) . Notice first that 9[^^). 


Define for shortness v'^ := §^,v'/ := obeys the recursion rule 


vl = -9,vl_, + - vt-1 + (24) 

obtained by differentiating the definition of vt- In particular is a linear function of 


vt-i,vt- 2 , • • • • This implies that vj-(9[^^} is uncorrelated with vt(9\^^) = 
holds for nf (@1^^). Using these properties, one can evaluate 




(«)) 


VE 






vtvl. 


= E 


, 2 




J_ _ V 

.2aj 2 u)'^. 



VV'E 






The matrix VV'E 


E 


d'^) 




E [-^v[vt] 


E 


1 

2lJ^ 


M _ 1 

CJ 


= 0 , 




The same 


(25) 


E 






is nonsingular unless vl-(9^^^) is zero a.s.; this cannot 


happen, otherwise from 


we would obtain a nontrivial relation among the at 


different t’s, but since is a white noise process with variance > 0 this is impossible. 


vm 


M\2 


exactly in 


2(TUT)2 


Indeed, with some transfer function machinery, one can evaluate E 
terms of the system parameters (see Lemma HD, but here it is enough to prove that it is 
nonzero. 

Notice moreover that the Hessian is positive definite; hence the negative log-likelihood 
has a local minimum, and the likelihood has a local maximum. 

Lemma 3. Consider the constrained maximum-likelihood estimator {9\^\(jj^'^'l) = argmaxL(0i,t<;) 
for a process x[^\ This estimator is asymptotically consistent. 
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Proof. This follows from standard maximum-likelihood theory, for instance from Ling & 


(w) 


has a local maximum in 


McAleer(2010), Theorem 1. We have proved above that E 

so using a suitable neighborhood of this point as the parameter space 0, the 
assumptions there are satisfied. □ 


Lemma 4. The estimator 


is asymptotically normal, i.e., \/iV(/3 — /3) —)• iV(0,H) in law. 

Proof. Once again we follow Poloni &: Sbrana(2014), Theorem 11. In view of their proof, 
it is enough to prove that the functions 


A 

ddi 


- E 


^ pi^) I ttO 


A 

doj 


- E 


^ fi'^) I -pO 


(26) 


(where the cr-fields T" are costructed starting on If = 
w G W. We recall that (cfr. ([25]) 1 


) decay with rate < 1 for each 


(27) 


d9i ’ dco 

The process vt{x^'^\9i) is a linear process in Yt. Using backshift operator formalism, we 
have 


Zt — Cf-2 Y ~ 2ef_i -|- ef_2 — [L — 2L -|- ^ 


= w[zt = w[\l — 2L + L^] 


6 




(.(»)..0 = 77(r)-‘T> = «(r)-.»;[7-2r + r= - 
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with H{L) = 1 — 6iL — This rational expression can be turned into a power series 

in L as in (1201) . with coefficients decaying as the maximum modulus of the roots of 
H{L), which are all smaller than 1 by the invertibility assumption. 

Similarly, one can evaluate 

with H'{L) = = —1 — Again, one turns this expression into a power series 

in L and obtains that the decay rate in the coefficients is exponential with rate given by 
the roots of H{L). Hence is also a linear function of Yj with exponentially decaying 
coefficients. Now we simply apply Lemma [T] to get the required decay properties for (1261) . 

□ 

Having proved Lemmas [3] and 01 we can conclude as in Poloni &: Sbrana(2014), Theo¬ 
rem 12 to prove Theorem [H 


H'jL) 

H{Ly 


[I-2L + L2 l2] 
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Figure 1: BELGIUM: Industrial production and associated smooth-trends 




The chart shows the industrial production series (in grey) together with the trends extracted with the 
standard HP filter (tiny black line), with the META approach (thick black line) and with the AMB 
approach (dotted line). The first chart refers to the sample 1974-1994, while the second one refers to the 
sample 1994-2014 
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Figure 2: FRANCE: Industrial production and associated smooth-trends 




The chart shows the industrial production series (in grey) together with the trends extracted with the 
standard HP filter (tiny black line), with the META approach (thick black line) and with the AMB 
approach (dotted line). The first chart refers to the sample 1974-1994, while the second one refers to the 
sample 1994-2014 
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Figure 3: GERMANY: Industrial production and associated smooth-trends 




The chart shows the industrial production series (in grey) together with the trends extracted with the 
standard HP filter (tiny black line), with the META approach (thick black line) and with the AMB 
approach (dotted line). The first chart refers to the sample 1974-1994, while the second one refers to the 
sample 1994-2014 
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Figure 4: ITALY: Industrial production and associated smooth-trends 




The chart shows the industrial production series (in grey) together with the trends extracted with the 
standard HP filter (tiny black line), with the META approach (thick black line) and with the AMB 
approach (dotted line). The first chart refers to the sample 1974-1994, while the second one refers to the 
sample 1994-2014 
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Figure 5: NETHERLANDS: Industrial production and associated smooth-trends 




The chart shows the industrial production series (in grey) together with the trends extracted with the 
standard HP filter (tiny black line), with the META approach (thick black line) and with the AMB 
approach (dotted line). The first chart refers to the sample 1974-1994, while the second one refers to the 
sample 1994-2014 


32 




Figure 6: PORTUGAL; Industrial production and associated smooth-trends 




The chart shows the industrial production series (in grey) together with the trends extracted with the 
standard HP filter (tiny black line), with the META approach (thick black line) and with the AMB 
approach (dotted line). The first chart refers to the sample 1974-1994, while the second one refers to the 
sample 1994-2014 
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Figure 7: SPAIN; Industrial production and associated smooth-trends 




The chart shows the industrial production series (in grey) together with the trends extracted with the 
standard HP filter (tiny black line), with the META approach (thick black line) and with the AMB 
approach (dotted line). The first chart refers to the sample 1974-1994, while the second one refers to the 
sample 1994-2014 
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Figure 8: UK; Industrial production and associated smooth-trends 




The chart shows the industrial production series (in grey) together with the trends extracted with the 
standard HP filter (tiny black line), with the META approach (thick black line) and with the AMB 
approach (dotted line). The first chart refers to the sample 1974-1994, while the second one refers to the 
sample 1994-2014 
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